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ABSTRACT 


An approximate method is presented which yields the supersonic 
nozzle profile necessary to keep an expanding high temperature gas mix- 
ture approximately in chemical equilibrium. From equilibrium calculations 
for given chamber conditions, the freeze area ratio is determined by use 
of a modified Bray criterion. A plot of area ratio versus a function of noz- 
zle geometry defines a limiting nozzle angle beyond which freezing is 
forecast to occur. This permit gan.exRag sion angle at each area ratio 
to be chosen to Keep the flow near equilibrium. The effectiveness of the 
resulting nozzle profile is evaluated by comparing exact calculations 
with equilibrium calculations. Contouring is shown to be effective for 
small nozzles when a correction factor is applied. The hydrogen-oxygen 


system is used as an example. 
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be, Introduction. 

During the last decade a marked increase of interest in rocket design 
and performance prediction has occurred. Payloads for space :exploration 
missions often differ greatly from those of a military nature, but the vehicles 
involved have a major common factor which is the rocket engine. Methods 
of propulsion currently under investigation for the future include ion, mag- 
netoplasma, and nuclear fission; however, the primary type today is the 
Chemical rocket engine. 

The idea of jet propulsion is old. In the chemical rocket, potential 
chemical energy of the fuel and oxidizer is converted into thermal energy in 
the combustion chamber. A convergent-divergent de Laval nozzle then accel- 
erates the gas stream from subsonic to supersonic velocity. This converts 
the thermal energy and pressure potential energy of the nearly stagnant, hot 
gas ahead of the nozzle into the kinetic energy of the supersonic stream 
which emerges with lower temperature and pressure. A reaction is imparted 
which is directly related to the time rate of change of the momentum of-the._. 
ejected matter. 

At the high temperatures engountered in modem combustion chambers, 

a large part of the chemical energy released may be used in breaking chemical 
bonds to produce simpler species (dissociation). As the temperature and 
pressure decrease during expansion, recombination of the simpler species 
occurs with a resulting release of chemically bound energy to the gas flow. 

A significant portion of the chemical energy is not available for conversion 


to kinetic energy unless reassociation of the dissociated combustion products 


| 


occur. Thus, the thrust produced by any propulsive system in which combus- 
tion products are expanded through a nozzle will be reduced if the recombinarc: 
tion process cannot occur rapidly enough to remain in equilibrium. 

In a rocket nozzle the elementary reactions between the dissociated 
combustion products seldom proceed at rates which are sufficient to main- 
tain true chemical equilibrium throughout the nozzle expansion fi 6]*. For 
a hypersonic ramjet, calculations indicate net thrusts which are attractively 
high if equilibrium is maintained, but which approach zero around mach num- 
bers of ten or eleven if the combustion products are frozen, i.e., constant 
composition [3]. 

To facilitate further discussion, the following descriptions of flows are 
given: 

Prozen: The chemical reaction rates are zero with composition 

Pee isnt Typical non-reacting flows come under this 
Classification. 

Equilibrium: The chemical reaction rates are infinite, and composition 
of the gas adjusts instantaneously to changes in temper- 
ature and pressure. 

Exact: The chemical reaction rates are finite. A finite time is 
required for the composition of the gas to adjust to the 
changes dictated by the pressure and temperature variations 


in the nozzle. 


* Numbers in brackets indicate references on page 30. 
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Near-equilibrium flow and near-frozen flow imply a small deviation 
from the artificial cases of zero or infinite reaction rates without an appreci- 
able change in the resulting performance prediction. Equilibrium flow yields 
maximum performance while frozen flow yields minimum performance. 


Consider the following typical three-body reaction: 


A+ B+M a AB + M 
R" 
where M is a third body which removes the energy during recombination, and 
R’ and R" are the rates of the reaction in the forward (recombination) and 
backward (dissociation) directions respectively. In the combustion chamber, 
the reaction would move to the left. At lower tempkratures and pressures, the 
equilibrium composition would shift toward the formation of AB. 

If, however, the reaction rates. are not fast enough to maintain the 
equilibrium conditions, recombination will lag. The gas leaving the nozzle 
will have some dissociated products which have not contributed to the kinetic 
energy of the flow. A reduction in thrust will occur if the recombination pro- 
cess does not occur rapidly enough to remain in equilibrium. 

It is assumed that the working fluid is a mixture of ideal gases; that 
the flow is steady, adiabatic, and quasi-one-dimensional; and that the effects 
of diffusion, heat conduction, and viscosity are negligible. Equilibrium flow 
calculations involve the simultaneous solution of the algebraic equations 
for conservation of mass, energy, momentum, and species plus the equation 
of state. Frozen flow calculations are the same except that the species 


composition is constant during the flow. These two types of calculations 
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present no particular difficulties. Many computer programs have been written 
to uyield gas dynamic and thermodynamic data at selected stations (4. SwGil i 

For an exact flow, though, the problem is more complex. Sarli, etal., 
[16] , summarizes the difficulties. The calculations require the simultaneous 
solution of the mass, momentum, and energy conservation equations; state 
and mixture enthalpy equations; kinetic differential equations for the finite 
rate of production of molecular species; and atomic continuity equations for 
all atoms present. The equations used in the exact solution are given in 
Appendix I. Again the assumptions are that the flow is one-dimensional, 
inviscid, and adiabatic, with.negligible diffusion. The simultaneous solu- 
tion of the above equations requires numerical techniques; hence, a digital 
computer is desirable [9, 10,16}. . 

Three basic problems exist, however, regarding these exact calculations 
[i 4. These are: 1. The numerical integration of the equations, when near- 
equilibrium conditions exist in the flow, is relatively slow and potentially 
inaccurate; 2. Itis necessary to specify the detailed mechanism and in- 
dividual rates for the reactions taking place; 3. The cost of the computer 
time involved prohibits casual investigations of complex, reacting flows. 

Today's space ventures require thrust prediction and/or weight minimi- 
zation such that a design based solely on a frozen flow calculation or an 
equilibrium flow calculation would be unacceptable. This does not mean 
that equilibrium and frozen calculations are to be abandoned. In view of the 
computational problems mentioned above that exist with the exact solutions, 


a method of establishing guidelines for the design of rocket nozzles is desir- 
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able which utilizes approximate methods. 

One of the first to investigate a modified equilibrium/frozen calcula- 
tion to predict actual performance was Penner [12]. Criteria were established 
which, in some instances, permitted a nozzle flow to be classified as near- 
equilibrium or near-frozen. Thus, fora given set of conditions, a reasonable 
prediction of performance could be given without encountering the complexi: - 
ties ofian exact calculation. Reichenbach [15] proposed an iterative proced- 
ure for the solution of nozzle flow problems. Agreement with exact calcula- 
tions was generally good. Several investigators da 11| have considered the 
design of an optimum nozzle. 

Another early investigator in this area was Bray [2. Si, His work has 
received considerable attention since its first presentation in 1958. Many 
of the current papers in the field of nozzle flow with chemical reactions 
utilize his concepts. 

The Bray criterion is an approximate procedure for prediciting where a 
reaction has departed significantly from equilibrium in a reacting nozzle flow. 
This prediction is made by determining the point at which the forward rate of 
reaction, R’, becomes of the same order as the rate required to maintain 
equilibrium. The flow upstream is considered in equilibrium, while that down- 
stream is considered frozen (invariant in chemical composition). While Bray's 
original work dealt with a single reaction system, later investigators (3, a 
8, 16] cextended the principle to include multiple reaction systems. 

An important point to consider is the effect of nozzle size and geometry 


on performance. Penner points out [12] that the reaction time must be small 
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“ompared with the residence time to achieve near-equilibrium. Hence, for 

small nozzles (short residence time) it is difficult to achieve near-equilibrium 
flows. 

After the above mentioned calculations are completed, it is then appa- 
rent whether or not the boundary configuration chosen (nozzle geometry) is 
likely to give the desired performance, i.e. thrust and specific impulse. A 

series of computations using different nozzle geometries would eventually 
yield a suitable nozzle design, but this would not be economical in terms of 
expense and labor. 

Kushida (8 ] and Koppang, et al., (7 presented a method utilizing the 
Bray criterion to determine the freeze point fora given nozzle. Reichenbach 
[1 4] suggested to this investigator that a reverse approach might be plausible. 
Since an estimate is available from Bray's criterion for the location of a freeze 
point as a function of area ratio, a nozzle designed to avoid the conditions 
that lead to freezing of the flow could approach the optimum performance for 
a given system. A further advantage would be the relatively simple calcula- 
tions involved which use equilibrium data. 

An analysis of this approach was undertaken in this paper. The system 
investigated was hydrogen and oxygen. To compare and evaluate the results, 
an exact calculation of the system was performed. 

This paper was prepared during the 1965-1966 academic year at the 
United States Naval Postgraduate School, Monterey, California. Acknowledge- 
ment is gratefully made to Dr. R.E. Reichenbach of the Aeronautics Department 


' 
forhis guidance and consultation as thesis advisor, and to the staff of the 
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computer facility for their cooperation and assistance in adapting the com- 
puter programs to the Control Data Corporation 1604 system. 
Bie Analysis. 

A common form g S, a6, 17] for Bray's sudden freezing criterion is 


the following semi-theoretical relaxation parameter: 


Equilibrium rate of change of species j (1) 


B = ‘Rate of the recombination reaction R; 

Thus, B relates equilibrium composition change to the exact composi- 
tion change. If the ratio is large, the reaction would be near-frozen; if the 
ratio is small, the reaction would be near-equilibrium. It is generally agreed 
that the freeze point is determined when B is near unity pe vj g|. 

For a three-body recombination reaction as shown in equation (1), the 
recombination rate, R', and the dissociation rate, R", are given by [7]: 

R' = kg n° Xp Xp Xyy | (2) 

R" = kp) n2X,,X), (3) 
where kg and kp are the reaction rate constants (assumed to be a function of 
temperature only); n is the molar density of the gaseous mixture; and Xn, Xp, 
Xp By and Xu are the mole fractions of the respective atoms and molecules. 

It has been shown[7, 8] that the bimolecular atom transfer, or meta- 
thetic reactions, are an order of magnitude faster than the three-body react- 
ions. Also, the energy released by the metathetic reactions is small com- 
pared with the three-body reactions. Consequently, the metathetic reactions 


are assumed to be in = equilibrium for this study. Equation (4) is a typical 


metathetic reaction. 
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R! 
OH+ 0 = H+ 0, (4) 
R" 


The main reactions to consider in the hydrogen-oxygen system are the 


following [7, 8, 16]: 


f 


l H+ H+ M= H, + M 


20+ 0+ M 0, + M (5) 


| 


2 


39 oH +30e+ @M OFlea» M 


40° OH + Het ir = H,O-+ M 
It has been found [7] that a system of this type is dependent on the 
dissappearance of hydrogen to determine the freeze point. As shown by .-. 
Kushida (8 ] and Koppang, etal. (7]. the net rate of disappearance of monatomic 
hydrogen is given by: 
dK, ,/m) 


-mn Fe — Ri + Rig + Rig eRe aeR's - By, (6) 


7 
where: mis the mean molecular weight, and the subscripts on R refer to 


the reactions in equation (5). 


The reaction rates are determined as shown by Penner [12] ; 


R', = ky (2)(nX,)(nX,)(nX\9 = 2 kynXX,, 
R'g = k3(1)(nX,) (nX)nX,) = kan Xe X Ky (7) 
Pe 


g = kgQ)(nXQ} (Xp (nXy) = kK gnSXOpXpXyy 


Since Xn is the mole fraction of any third body, itis unity. Rewriting 
equation (7): 


R! = 2 keen ge (8) 
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To obtain the expression for B as shown in equation (1), the equilibrium 
rate of change of hydrogen concentration is divided by the rate of the recom- 


bination reactions for hydrogen. It should be noted that R' represents the 


rate of recombination and R" represents the rate of dissociation. 


3 = mnd(H/m)/dt 


R') + Rig t+ R'y 


_ mnd @H/m)/at 


2, n9Xpq? + kgn3XpXo + kynSXouXy 


_ mdH/m)/dt (9) 


Xpqn2(2kyXqy + k3XoO + kgXoH) 


Be. d{ln m dt 
n2(2 kyXpy + kgXq + kyXon) 


Introducing v = dL/dt and dividing both sides by d(@) where € =A/A* 


(area divided by the throat area): 


B _ vgin XH/mjl v/a 


de n@(2 kyXp + kX + k4yXop) de 
Rearranging: 
B _ vdfin XH/m)] fae 18) 


d@/db 2 (2 ky Xqz + kgXG + kgXoOq) 

The right side of equation (10) can be calculated using equilibrium 
data fora given system. The rate constants, ke and kj, are assumed to be 
functions of temperature only. By plotting the right side of equation (10) 
versus area ratio, a means of determining the freeze point is provided. For 
every freeze area ratio, there is a corresponding value for B/(dé/dL). 


(Curve A, Fig. 1) 






Equilibrium 
Margin 
(EM) 


dé /dL 


Freeze Area Ratio 


Figure 1: Freeze Point Location 
Since B is assumed to be constant and equal to unity, it is thus seen 
that the local value of dé/dL determines the freeze point. Therefore, if the 
nozzle is designed such that d€é/dL is less than the local value for each 


freeze area ratio (Curve A’), the flow should remain in near-equilibrium. If 





the nozzle is designed such that dé/dL is greater than the local value for each 


freeze area ratio (Curve A"), the flow should become near-frozen. 

The vertical distance between Curve A and Curve A' is arbitrary and 
not necessarily constant. It is important to note, though, that a small dis- 
placement from Curve A possibly would not yield near-equilibrium flow be- 
cause of the inherent inaccuracies introduced by the assumptions. A large 
displacement would increase the nozzle length because of the slower expan- 
sion to a given area ratio. This displacement will be referred to as the 
equilibrium margin (EM). Since the curves in Fig. 1 are usually plotted on 


semi-logarithmic paper, the constant displacement EM becomes a multipler 


to the basic curve. 


The information available from Fig 1 thus provides a guide for nozzle 


design which incorporates the effects of chemistry. 
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3. Procedure. 

The first step necessary in determining a nozzle profile is the calcular. 
tion of the right side of equation (10). All of the values used are those ob- 
tained from equilibrium flow with the exception of the rate constants. The 
rate constants, however, may be assumed to be functions of the equilibrium 
temperature at each point in the nozzle. Since equilibrium flow is indepen- 
dent of the nozzle geometry, the same values calculated for one system 
(chamber conditions, O/F ratio, and reactants) can be used in the study of 
any nozzle geometry, regardless of size. 

Values for d [In (p/m) | /dé may be obtained graphically or by analyt- 
ical. methods onge Xz, and m are known for different area ratios. This per- 
mits the right side of equation (10) to be calculated and plotted versus freeze 
area ratio, and now represents B/(d&é /dL) as a function of freeze area ratio. 

To convert this into an axisymmetric nozzle profile, the following 


geometric relationship may be used [7] : 


dé _ d(r/r*)* _ 2r d(r/r*) (11) 
dL dL rx dL 
and dr/dL = tan 8 (12) 


where 6 is the half-angle of the nozzle. Then 


dé _ 2r tand® . 2Vé tane (13) 
dL re op* r* 

ortang = dé —ch 14 
dL 2Ve sa) 


Since dé/dL is available from the previous calculation of B/(dé /dL), 


equation (14) provides local values of 6 for every freeze area ratio. It should 
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be noted that nozzle size is now present~ in that equation (14) involves the 
throat radius, r*. 

The half-angle, 98, is conveniently calculated by obtaining an expres- 
sion for B/(dé/dL) as a function of area ratio. If B is assumed to be unity, 
then a polynomial form for dé/dL is 


e4 (15) 


= ae, 3 
dL/dé = bh, + bj¢ + b,€ je bo a by 


where b,, bi. etc. are appropriate coefficients. 
Equation (15) may be integrated with the boundary conditions that 


G&= 10atL = OO. 


Z 3 4 fs) 
b E b,€  b,é 
L=b€E + ra aa en BoSen. + 0 a (16) 
O Z 5 4 3 


Using equations (14), (15), and (16), the half-angle and nozzle length 
may be calculated for each oe area ratio. Since the throat radius has 
been set, the throat area is known. With the area ratio and the throat area, 
the local area and thus the nozzle radius are determined. If an equilibrium 


margin is desired, it must be applied to both equations (15) and (16). 


Special provisions must be made near the nozzle throat. The right side 
of equation (10) becomes very small as the area ratio approaches one. Thus, 
dé/dL is large and, as equation (14) shows, the half-angle approaches 90 
degrees. This is unacceptable from a gas dynamics approach. While this 
difficulty is not the concer of this study, it should be recognized as an 
area requiring special handling. 

Some method of fairing in the profile near the throat must be used until 
the area ratio permits the contour determined by equation (14). When this 


point in the nozzle is reached, the values of L must be adjusted since they 


ZZ 





no longer conform to the calculated contour. This correction is a constant 
shift of the L values. As long as the faired section near the throat has local 
half-angles less than that given by equation (14), the flow should remain in 
near-equilibrium. Since the allowable half-angles in this region are large, 
no particular difficulty is encountered. 

At large area ratios, 9 becomes very small. Thus, it is usually desir- 
able to terminate the contouring after the majority of the benefits of increased 
recombinationihnave been realized. Further expansion may be accomplished 
at more rapid rates with consequent near-frozen flow. 

4. Results and Discussion. 

A study of a particular system was performed to check the validity of 
the proposed method of nozzle design. The hydrogen-oxygen system was 
chosen because of its relative simplicity and limited number of combustion 
products. A chamber pressure of 300 psia and an oxidizer-to-fuel weight 
ratio of six was chosen as representative of current combustion conditions. 

Equilibrium and frozen values for the above mentioned conditions were 
computed with the aid of a computer program written by Gordon and Zelznik 
& 6 | . Tabulated results may be found in Appendix II. The computer pro- 
gram was written for the IBM 7090 system. Since the computer at the Naval 
Postgraduate School is a CDC 1604, several modifications due to system 
differences were necessary so that the equilibrium calculations could be 
made. The basic program was unchanged. 

A least-squares curve fitting routine was used to obtain a third order 


polynomial relating In (X;,/m) and area ratio. 
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in (X,,/m) = ag tat a,é? + eee (17) 


where: aoe= 522623 


5.22515 x 1072 


aT = 
a, = -4.27257 x 107 
a, = 1.61956 x 107° 
Then d [in (X}/m)| /dé = a, + 2a7€ + 3a,€? (18) 


The rate constants kj, Ky, and ky were calculated from the following 


relationship: 
k = DT exp (F/T) (19) 
where T is the local equilibrium temperature and D, 5, and F are constants 


listed in Table I [18]. 


TABLE I 
REACTION RATE CONSTANTS 


Reaction D E F 


l 2.3555x 1017 -1.4998 -1493 
3 5.5000 x 1o!! 0.0000  -10000 
4 1.7268x 1045  -1.0900 5375 


Using equation (18), equation (19), and the equilibrium values, the 
right side of equation (10) was then calculated. The results are given in 
Table II. As seen from equation (10), these values also represent B/(dé/dL). 
Thus, the lecus of points shown in Fig. 2 (EM = 1.0) defines the hypothesized 
boundary below which departures from near-equilibrium flow probably occur. 
The values for EM = 1.5 were obtained by multiplying the basic curve 
(EM = 1.0) by 1.5. The choice of 1.5 was arbitrary. In view of the 


possible inaccuracies introduced by the simplifying assumptions made in 
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the formulation of the problem, it was assumed a factor of approximately 
EM = 1.5 was sufficient. It should be remembered that EM = 1.0 represents 
the point where, theoretically, the flow is just able to maintain near-equili- 
brium. A_large value for EM would result in a long nozzle since the expan- 
sion would be slower. Also, an additional important aspect of this design 
method is to determine the maximum rate of expansion which may be used 
and still maintain near-equilibrium conditions. 
A polynomial of B/(d € /dL) as a function of area ratio was obtained 
with a curve-fit routine which resulted in the following coefficients: 
b, = 9.28458 x1073 
b) =-2.39544 xl0~ 
by = 1.34259 x1074 
b, = 3.95404 x107° 
b, =-3.83076 xl0™° 
Equation (14) was used to determine the maximum expansion angle 
that the flow could tolerate without departing from near-equilibrium. These 
values are shown in Table III. It should be noted that small area ratios yield 
large values for @. 
Calculations were performed with nozzle throat radii of 2.0 inches, 
0.96 inches, and 0.12 inches. These different size nozzles were identified 
as Run 1, Run 2, and Run 3. A designation was added to each run to indicate 
whether the nozzle was a straight conical expansion (S), a contoured expan- 


sion (C), or a contour expansion with EM = 1.5 (CM). For example, the run 


with the two-inch throat that had a contoured expansion with EM = 1.0 was 
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designated Run I1-C. 

A maximum of 45 degrees was used for all runs except Run 3=CM where 
8 was limited to 30 degrees. In all cases, a circular arc of radius 0. 4 times 
the throat radius was used to fair in the profile near the throat [13]. A Rin- 
ear section was used on the contoured nozzles (1-C, 1-CM, 2-C, and 
3-CM) which was carried to the point where the area ratio indicated that the 
contoured expansion angle as shown in Table III was permissable. The 
straight conical expansions were all at 15 degrees. See Table IV fora fur- 
ther description of the nozzles. Figs. 3, 4, 5, and 6 show the profiles used. 
Tables V to XI give the nozzle geometry for each run. 

Exact calculations for all seven runs were made with the computer 
program written by Zupnik, Nilson, and Sarli fig}. As with the equilibrium 
computer program, certain system modifications were necessary to convert 
from the IBM 7090 to the CDC 1604. Tabulated values obtained in these cal- 
culations are shown in Appendix III. 

The benefits of maintaining equilibrium flow can be seen from Fig. 7. 
Even at the throat, the specific impulse difference is five secs. As the ex- 
pansion continues, the difference increases to 20 secs at an area ratio of 40. 

Figs. 8, 9, and 10 illustrate the effect of contouring on performance 
for the different nozzle sizes. The contour is most effective on nin 3-CM 
and is insignificant on run 2-C. The apparent explanation for run 2-C is 
that this run is not modified. That is, it does not have the equilibrium margin 
that runs 1-CM and 3-CM have. Also, a small nozzle is subject to earlier 


freezing than is a large one. This is due to the smaller residence time ina 
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small nozzle which gives less time.ior the recombination reactions to take 
place. 

Runs 2-C and 2-S illustrate the effect of any contour on nozzle length. 
Fig. 12 shows that a specific impulse of 365 lb-sec/lb may be obtained with 
a 0.10 foot nozzle for the contoured nozzle. The straight conical nozzle, 
though, needs 0.17 feet for the same performance. Fora given langth of 
0.08 feet, the specific impulse difference is 28 secs. 

The most significant result is that shown by the small nozzle. Even 
at relatively small area ratios the effect of contouring is evident. At an area 
ratio of 40, the difference in specific impulse is six secs. While this is not 
a large percentage increase, the effect over a long period would be appreci- 
able. Thus, long burn times would require less fuel for the same total per- 
formance. Where weight is a critical factor, contouring is particularly ad- 
vantageous. Fora hotter system such as hydrogen=fluorine, the percentage 
increase in performance is expected to be more significant. 

The variation of hydrogen concentration is given in Figs. 13, 14, and 
15. The equilibrium value is seen to decrease uniformly with area ratio. 

In each nozzle size, though, the exact concentrations show a leveling off 
which indicates a departure from equilibrium. For the conical nozzles the 
change is gradual, but it occurs where the allowable expansion angle is 
exceeded (15 degrees). All of the contoured nozzles showi: a marked "knee" 

in the hydrogen concentration when the contouring is stopped (aboutan area:ratio 
of 4.0). Thus, the effect of contouring a nozzle to delay departure from 


equilibrium appears to be significant. 
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Each of the figures illustrating hydrogen concentration shaws.a gradual 
deviation from the equilibrium values at small area ratios. It is felt that this 
is the effect of finite rate constants and not due to inaccuracies in the exact 
calculations. When large rate constants (k = 1029) were used, the exact 
values essentially duplicated those of the equilibrium calculations. Also, 
when small rate constants were used (k = 100), the concentrations were 
invariant indicating frozen or non-reacting flow. 

The temperature distribution during expansion is given in Fig. 16. 

Run 3-CM has larger temperatures than Run 3-S. This indicates that more 
recombination is occurring in the contoured nozzle than in the straight con- 
ical nozzle. 

As mentioned above, the amount of equilibrium margin used was arbi- 
trarily chosen as 1.5. Further investigation should indicate an optimum 
factor for EM. A further improvement would be to correlate the increase in 
weight due to a longer nozzle with the increased performance. 

The approach used in this thesis was formulated with hydrogen dis- 
appearance as the indication of the freeze point and subsequent deviation from 
near-equilibrium conditions. Further study might show that other species 
such as OH or monatomic oxygen give better results. Possibly a combina- 
tion of species could be used. Sarli, etal. [16] investigated this modified 
approach in determining the freeze point for multiple reactions. This could 
serve aS a significant refinement to the basic method of determining optimum 


nozzle contours. 
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De Conclusions. 

Optimum contouring of a nozzle by use of a sudden freeze criterion 
results in a delay of the departure from near-equilibrium flow and an increase 
in performance. This is relatively more effective with smaller nozzles where 
recombination normally lags because of the short residence time. 

Direct application of the Bray criterion for sudden freezing does not 
yield the best results; however, an equilibrium margin applied to the basic 
curve increases the performance. Further study in this area should indicate 
an optimum equilibrium margin. 

Hydrogen may not be the ideal species to indicate the freeze point 
in determining the optimum nozzle profile for this system. Other possibilities, 


including combinations of species, should be investigated. 
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. 083 
124 
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32009 
353 
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TABLE II 


FREEZE POINT LOCATION 


B Area 

qd /dL Ratio 
0036307 1.903 
.0045955 2.079 
. 0056646 PANY | 
. 0068524 2.420 
. 0081592 3.976 
.0111240 5.902 
- 0184950 7.891 
. 0277860 13eS00 

TABLE Ill 


MAXIMUM EXPANSION ANGLE 


Half-angle (Degrees) 


Run Run 
1-CM 2-C 
86.54 85.20 
84.01 SO 
74.79 G9. ol 
53. 86 44.59 
44,89 B32 69 
33. 46 29.45 
30. 40 22.90 
Zone 18.66 
20.58 16.13 
14.13 Os 27. 
9=02 6.5/2 
4.28 3.08 
2.05 1.47 


B 
d_/dLb 





- 039031 
052248 
067470 
. 084753 
- 365240 
1.102000 
2.687200 
39. 807000 


Oooo Ff kK ND PD W 
INO 
WO 


Run 


1-S 


1-C 


(Fon ® 


3-5 


3-CM 


Throat 
Radius 


(ft) 


0. 167 


0.167 


OPE G/ 


0.08 


0.08 


0.01 


0.01 


TABLE IV 


NOZZLE PARAMETERS 


Throat 
Area 


_ét2) 


0.0873 


0.0873 


0.0873 


0.0201 


0.0201 


0.003142 


0.003142 


Circular arc to 15° straight expan- 
sion. 


Circular arc to 459 ; then 45° straight 
to 45° allowable contour; then con- 
toured with EM = 1.0 to point of 
coincidence with Run 1-S; then 159 
straight expansion. 


Same as Run 1-C except contour is 
with EM = 1.5. 


Same as Run 1-S except different 
size throat. 


Same as Run 1-C except different 
size throat. 


Same as Run 1-S except different 


size throat. 


Circular arc to 30° ; then 30° straight 
to 30° allowable contour; then con- 
toured with EM = 1.5 to area ratio 
of 4.0; then 15° straight expansion. 


Radius 


ee 


Oe 6/ 


O21 9 


0.226 


0.240 


0.253 


Oe Ze 1 


0.288 


0.300 


OR, 


0.333 


0.351 


1.042 


TABLE V 
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NOZZLE GEOMETRY FOR RUN 1-S 


0.3484 
0. 3879 
0. 4631 
0. 668 

0.8254 
1. 2634 
1.7438 
2.2734 


3.4081 


Zoe 


39. 


OZ 


» 444 


. 306 


az 


458 


» 476 


~982 


050 


U5Z 


0.328 


0.685 


0.840 


1.002 


1. 198 


1.563 


2.276 


2.727 


3.692 


Radius 
ip) 


0.167 
GO. ZUS 
0.234 
0.249 
0.267 
0. 283 
0.301 
0. 309 
0.322 
0.334 
0.354 
0.395 
0. 439 
0.492 
0.589 
0.780 
0.901 


1.160 
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TABLE VI 


NOZZLE GEOMETRY FOR RUN 1-C 


Area 


(ft2) 


0. 0873 


0.1504 


0.1715 


0.1950 


0.2234 


0.2514 


0. 2839 


0.3003 


0.3251 


0.3511 


0. 3943 


0. 4896 


0.6052 


0.7606 


1. 0907 


1. O32 


2.9914 


4.2247 


Area 
Ratio 


1. 000 
1.723 
1,965 
2.234 
2.960 
2.881 
3.253 
3. 441 
3.726 
4.023 
4.518 
9.610 
6.934 
8.715 
12. 497 
21.923 
29.235 


48. 410 


Length 


0.000 
0. 080 
0.093 
Oo EOZ 
0.115 
On ts7 
0.182 
0.312 
0. 449 
0.640 
0.701 
0.837 
Ue Site 
Us ales: 
1.836 
2.699 


3.991 


Radius 


36 


TABLE VII 


NOZZLE GEOMETRY FOR RUN 1-CM 


Area 


(ft2) 


0.0873 


0.1507 


0.1694 


0.1804 


0.1936 


0.2108 


0.2366 


0.2838 


0.3169 


0.3539 


0. 3897 


0. 4740 


0.5609 


0.8005 


We ee Wy 


2.4886 


4.0184 


Area 
Ratio 


1.000 
weal 
1.941 
2-106 
2.200 
2.414 
Poe (El yl 
3.251 
3.631 
4.055 
4.466 
5. 431 
6, 427 
3, Les 


15.594 


28.916 


46. 046 


TABLE VIII 


NOZZLE GEOMETRY FOR RUN 2-S 


Length Radius Area Area 

i? ie” (ft2) Ratio. 
0.000 0.080 0.0201 1.000 
0.023 0.085 0.0226 1.124 
0.035 0.088 0.0245 Ls 27 
0. 069 0.097 0.0298 1.481 
0.103 0.107 0.0356 l 72 
0.149 0.119 0.0442 2.198 
0.170 0.124 0. 0486 2.417 
0.198 0. 132 0.0546 2.717 
0.246 0.145 0. 0658 3.270 
0.275 0.153 0.0731 3.636 
0.312 0.163 0. 0826 4.125 
0.349 0.172 0.0932 4.635 
0.422 0.192 0.1155 9.746 
0.505 0.214 0.1440 7.160 
0.673 OF Zi59 0.2108 10. 482 
0.916 0.324 0. 3303 16.422 
1,255 0.415 0.5412 26.911 
1. 463 0.471 0.6963 34.625 
1.631 0.516 0.8356 41.552 


37 


TABLE IX 


NOZZLE GEOMETRY FOR RUN 2-C 


Length Radius Area Area 

aie) ef) cree Ratio 
0.000 0.080 0.0201 1.000 
0.017 0.085 0.0226 Pe Acc) 
0.024 0.03% 0.0260 L.295 
0.030 0.097 0.0296 1.473 
0.038 0.105 0.0346 L722 
0.051 0.114 0.0414 2. 046 
0.067 OalZZ 0.0465 2.314 
0.095 0.130 0.0527 Zoe. 
0.156 0.140 0.0616 3.062 
0.234 0. 149 0.0693 3.445 
0.308 0.161 0.0819 4.073 
0.356 0.175 0.0958 4.764 
0.442 0.198 O.1 225 6.093 
0.624 0.246 0.1895 9.425 
0.801 0.293 0.2700 13.424 
0.981 0.342 0.3667 18.236 
Las2e 0.432 0.5875 29 5200 
Looe 0.488 0.7486 37.160 
1.669 0.526 0.8689 43.092 
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Length 
x 102 


(ft) 
0. 000 


0.204 
0.432 
0.808 
1. 607 
2.368 
2-992 
3.935 
9.991 
9.933 
14.211 


18. 490 


Radius 
x 1024 
(ft) 


1.000 
1,041 
1.101 
17202 
1.417 
1.620 
1.787 
2.040 
2. a1 
3.647 
4.794 


9.940 


TABLE X 


NOZZLE GEOMETRY FOR RUN 3-S 


Area 
x 103 


ti) 
0.3142 
0.3402 
0.3810 
0.4538 
0.6308 
0.8246 
1.0035 
1. 3079 


2.1093 


4.1794 


7.2194 


11.0843 
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TABLE Xl 


NOZZLE GEOMETRY FOR RUN 3-CM 





Length Radius Area Area 
x 102 x 102 x 10° Ratio 
(ft) (ft) _(ft2) 
0.000 1.000 0.3142 1.000 
0.176 1,041 0.3402 1.090 
0.371 els 0. 4023 1.28! 
0.563 Ney dete. 0.4400 1.401 
Lad 1.334 0.5593 1.780 
3.199 1,423 0.6358 2.024 
6.940 lias 0.7576 2.412 
12.5595 1.668 0.8738 2.781 
29.700 Lagug 1.0289 Cy AIS 
CE). GIS ae. 1.1547 3.676 
52.417 2.006 1.2639 4.023 
2/ UZm 2.147 1. 4488 4.612 
97.671 Zee) 1.6981 5. 406 
59 9136 2.985 2.7995 aig eli 
63.627 3.961 4.9299 Soa 
Do. ooo 4.705 6.9552 Mido Meas, 
68.324 Sea 2 8.5669 Zi 20 
70.170 5. / 17 10. 2664 CYA OS, 
720s 6.211 LZ. Lio 38.576 
72.936 6. 458 3030 41.708 
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APPENDIX I 


EQUATIONS FOR EXACT CALC ULATIONS 


The equations involved in the analysis of nozzle flow with a chemical 
reaction include those of gas dynamics and of reaction kinetics [18] : 
Gas Dynamic Equations 


Overall mass conservation: 


fAw = mn" (I-1) 
Total energy conservation: 
S ee 
= hi Te +t igT = Hp (I~2) 


where: 


- 
hs “f Cr: dT + (He: ), 
Momentum conservation: 
ada + GedP = 0 (1-3) 
Equation of state: 
P=PRTEAT (I-4) 


Speed of sound: 


| 
C= PTC) gh Ie (I=5) 


h 
= [gvrT2u% ]* a6) 
where ¥ is defined as: 
ee de _ (I-7) 
Mach number: 
M= 4% (I-8) 
C 
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Reaction Kinetics Equations 
General form of im. chemical reactions: 


oy. Ree “ y- 
4 @ e Silex, ® 2 & 
Sj Gy “ (1-9) 


where ZA; and Zz are the stoichiometric coefficients of the jth species 
in the jth reaction, reactants and products respectively, and K gj and Kpj 
are the specific reaction rates of the forward and backward chemical changes 
in the jth reaction of the chemical system. 

The rates of change of the concentration of the species M, are deter- 


mined by the relation: 


. Mr ” U ad 
EMT) = 2 (M4 - 4) {he I CM) 


(I=10) 

VK 

-R,, ir CM] } 
Ka | 
The atomic species continuity equations are: 
n 

Ge. s &:Ty = constant 

<eTr s2Pet | 44°39 (I-11) 


(4 = 1,2) +616y M9) 


where Bij is the number of the it® atomic species in the jth molecular species, 
n designates the total number of chemical species, and n, is the number of 
atomic species. 
The molecular species continuity equations are: 
n ‘, 
pool = B(x - 24) [hy BPR) | 

, ag (I-12) 
— Rey TT (PR) | ae 


ti 


ey 


It should be noted that of the n equations represented by (I-11) and 
(I-12), only n-l are independent because of the overall mass conservation 
equation (I-1). The exact calculations involve the simultaneous solution 


of equations (I=-1), (I-2), (I-3), (I-4), (I-10), (I-11), and (I-12). 


28 


APPENDIX II 


EQUILIBRIUM AND FRQZEN RESULTS 


TABLE II-1 
Equilibrium Composition Values 


39: 


Area Pressure Pressure Temp Gamma Specific Enthalpy 
Ratio (Aymos) Ratio Deg K Heat Cal/gm 
(P./P) Cal/gm=-°K 
Chamber 20.414 1.000 3439 L.23f3 2.9563 0 
Throat 11.811 1.728 3274 1.1294 2.7698 -277.8 
1.083 8? 165 2.500 3166 1.1288 2.6193 -455.7 
1.124 7.423 227,50 3139 1.1288 2057 #3 -500.4 
] 268 6. 805 3.000 3114 1.1288 2.5380 -540.8 
Las 6.281 $. 200 3091 1, 1288 2.9010 -977.6 
1, 269 9. 832 3.500 3070 1.1289 2.4660 =611.3 
lass 5.103 4.000 3033 L. 1291 2.4013 -671.5 
Lb, 589 4.083 5. 000 2971 1.1297 2.2887 -770.0 
1.723 3.402 6.000 2920 1.1304 2iS'S2 -848.8 
1... 99'S 2.916 7.000 2877 L.tsr3 Zee 0S -913.8 
2 .iDiZ9 2. §52 8.000 2840 1.1322 2:00:37 6 -969.3 
2 i 2.268 9.000 2807 1. d632 Lagi7 27 017. 
2.420 2.041 10.000 2778 1.1342 1.9144 -1060 
3.976 1.021 20.000 2582 1.1442 1.5364 -1327 
a. 02 0.600 34.020 2425 1. 1567 1.2786 —Voul7 
7.891 0.408 50. 000 2306 1. 1682 ees -1646. 
13.399 0.204 100.000 2083 Leis i 3 Org 209 -1861 
45.921 0.041 500.000 1973 Levan. 1 0.7519 =2279 
70 009 0.020 5000. 000 989 1.2437 0.7190 -2419 


TABLE II-2 


Equilibrium Composition Values 





Area Mean Mach Ore Cp Specific Impulse 
Ratio Mol. wt. Number  Ft/sec (lb-sec/1b) 
Isp)y Isp 


ome SS SR ee Et Sa ee ee ee ee eee eee 


Chamber 13.043 0.000 Pe a ey = 


Throat Lo22.3 1.000 7659 ~654 29 3 159.95 
1.083 13.341 1. 308 7699 . 837 302.2 liom 
1.124 13. 370 M8 7 7655 . 877 305% 9 208.7 
1.168 NS es 8 |) 7g 1.440 76595 2 Z 309.6 216.9 
1,213 13.421 1.495 76595 ~942 313.0 224.2 
1.299 13.444 1,545 7655 969 316.3 230.6 
Liedio3 13.484 a seme 76959 1.016 3 2iiZ. 241.7 
1.539 bp. 548 1.769 76995 1.088 3373) 298.8 
Mwai 23 Lomo9.9 1.876 76595 1.142 340.0 2 Tlgeaf 
1.903 13.641 1.963 76595 1.1895 346.6 282.0 
Zn07 9 13.677 2.037 7659 lea! Jez. Z 290.4 
2.291 oer 07 2.101 76595 lao 357.1 297.6 
2.420 13.733 Paes yt 76595 eas. SOaro 303.7 
3.976 13.887 2.906 7655 1.429 Che. 339.9 
9. 9Q2 Ls 7 8 2.799 76995 AL 404.6 363.3 
7.409). 14.028 Nye eh One: 7695 Layo) 416.1 378.5 

13.355 14. 082 3.261 7655 L632 434.2 402.95 

45.921 14.111 4.085 76595 1.870 466.6 444.9 

77.299 14.112 4.479 76595 Leo28 477.2 458.8 


TABLE II-3 


Frozen Composition Values 


Area Pressure Pressure Temp Gamma 
Ratio (Atmos) Ratio Deg K 
ee ee ee! 
Chamber 20. 414 1.000 3439 1.2057 
Throat 11. 483 law78 31is5 1.2096 
1.065 8.165 2.900 2936 1,2i23 
1.101 7.423 2.750 2887 1.2130 
1.139 6.805 3.000 2843 1.2137 
1.179 6.281 3.250 2803 1.2144 
1,219 9. 832 Sm508) 276// 1. 2s 
1.301 S03 4.000 2702 1.2162 
1.464 4.083 9.000 2597 1.2183 
1.624 3. 402 6.000 2513 1.2201 
Lan78 2.916 7.000 2444 1.2216 
1293) Aes) Sh 8.000 2385 1.2230 
2.078 2.268 9.000 2334 1.2243 
2 aee2 2.041 10.000 2290 1. 22555 
3.523 NGO 20.000 2012 1.2329 
9.10918 . 600 34.023 1817 1.2413 
6.702 - 408 90.000 1684 1.2472 
11,048 . 204 100.000 1464 1.2590 
35. 866 041 500. 00 1035 1.2924 
99.600 020 1000.000 #£882 1.3074 


Specific 


“Heat 


m-OK 


Enthalpy 


— Cal/gm 


TABLE II-4 


Frozen Composition Values 


Area Mean Mach C* Cr Specific Impulse 
Ratio Mol. wt. Number  Ft/sec (lb-sec/1b) 
ee .LLLLL!l Isp)v Isp 
Chamber 13.043 0 ==55 a= = oe paar ote 
Throat 13.043 1.000 7 473 . 680 28057 138%" 
1.065 13.043 1.280 7473 . 846 25510 136. 6 
Labor 13.043 1.390 7473 . 886 298.8 ZU 
Las 13.043 1.412 7473 220 S110) aie) 213.6 
Ne TAS) 13.043 1. 468 7473 . 949 304.8 IATA Shoe 
Ne ike 13.043 lego re 7 473 .976 307.6 226.6 
1,301 13.043 1.606 7473 1.021 oZio Za7. 
1.464 13.043 1.748 7473 1.090 SZ 753. 1 
1.624 13.043 ses, 7473 1.141 RAL cus, 265.0 
12779 13.043 1 S/S) 7473 L eel 333.4 274.4 
i eu 13.043 PAR OCA 3: 7473 1.214 Scout 282. ! 
Zar 13.043 Use 7473 Ye 242 342.2 288.5 
Larne 2 13.043 Zino 0 7 473 1.266 345.7 294.1 
S2Zz0 13.043 2.941 7473 1.403 366.8 325.9 
9.098 13.043 2.830 7 473 1. 489 380.8 346.0 
6.702 13.043 oie Us 7473 1,544 SCG. 3 358.6 
11.048 13.043 3.422 7473 1.628 403.9 378.2 
35. 866 13.043 4.373 7473 L773 428.4 411.8 
99.600 13.043 4.830 7473 Tous 436.1 422.3 
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TABLE II-5 
EQUILIBRIUM FLOW SPECIES CONCENTRATIONS 


Mol Fractions xl 03 


Area QO} Hj Oo H, OH H,O 
ae! ll OU ee ee eee? Le Ce ... F_ 
Chamber 8.16 56.88 7.75 252.15 62.82 612.24 
Throat 6.06 49.35 6.14 248.38 52.17 637.90 

1.083 4.81 44,33 5.07 246.06 45.20 654.53 

1.124 4.51 43.03 4.80 245.50 43.44 658.71 

1.168 4,25 41.86 4.56 245.00 41.84 662.49 

1.213 4.01 40.77 4.34 244.56 40.38 665.93 

1.259 3.80 39.78 4.14 244.16 39.04 669.08 

1.353 3.43 37.98 3.79 243. 48 36.66 674.67 

1.539 2.85 34.99 3.22 242.44 32.76 683.73 

1.723 2. 43 32.56 2.79 241.69 29.67 690.87 

1.903 2. 10 30.52 2.44 241.13 27.12 696.69 

2.079 1.83 28.76 2.45 240.72 24.98 701.56 

2.161 1.61 27.22 1,92 240. 40 23.44 705/72 

2.420 1.43 25.85 1.72 240.15 21.53 709.32 

3.976 0.55 17.19 0.69 239.55 12.16 729.86 

5.902 i 4 11.26 0.26 240.17 6.74 741.36 

7.891 0.08 7.63 0.11 240.99 3.94 747.24 
13.355 0.01 -3.04 0.01 242.56 lel? 7G 
45.521 ---- 0.10 -a-- 243.94 0.01 755.95 
77.299 ---- 0.01 ---- 243.99 ---- 755.99 
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APPENDIX III 
EXACT RESULTS 
TABLE III-] 
EXACT FLOW SPECIES CONCENTRATIONS FOR RUN 1-S 


Species Concentration 
(lb-mol/Ab x 104 


Area O H O, H, OH HO 
a a ee 2! ee ee ae 
anaes 1.79 23.91 2.05 177.83 21.83 507.98 
1.834 Lae 23.66 2 902 177s 21.19 500875 
2.068 ] weil 22.45 er 176.51 19.27 511.30 
2.297 1933 dhe SY 1.63 175.72 17.68 513.43 
2.645 ya 20.10 1.41 1 7ART 7 1507 515.96 
2.979 1.00 now) 1.26 L7AmO7 14458  SD7MS2 
39037 0.92 18.69 lags 173.63 13.46 518.96 
3.662 0. 82 18.07 1. 06 173.03 pA 27) 520. 48 
3.992 OF77 17.73 1.00 172.64 11.53 5 Ze 1 
4,444 0.71 74) 0.93 72's 10.69 52245 
5.306 0.63 70 0.84 177 9.48 523.91 
7 ane 0.52 7s 0872 170.10 7.48 S2ene5 
9.458 0.47 17. 46 0.75 169. 33 6.52 527.35 
14.476 0.40 18. 42 0.61 167.81 4.86 529.23 
19.982 035 19.29 0557 16se70 3877 530. 45 
26.050 0.31 20.01 0.54 165. 84 2.99 531.33 
39, 052 0.24 203 0.51 164.70 1.98 5oam47 
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TABLE III-2 
EXACT FLOW SPECIES CONCENTRATIONS FOR RUN 1-C 


Species Concentration 
(lb-mol/lb) x 104 


Area O H O H OH H,O 
Te ee a. — - a ~~ 
1.723 1.79 23.91 2.05 177.83 21.80 507.98 
1.965 a7 23.93 2.03 177.06 20.57 509.37 
2.234 1.64 24.00 1.99 176.21 19.25 510.84 
2.560 1.44 22.90 1.81 e.g 86174 518,23 
2.881 1, 22 21.35 1.55 174.46 15.66 515.71 
3.253 1.00 19.51 1.28 173.69 13.84 518.31 
3.44) 0.90 18.66 1.16 178.37 13.03 529,46 
3.726 0.78 17. 49 1.00 172.97 11°99 520,99 
4.023 0.67 16.43 0.87 172.64 10.93 522.36 
4.518 0.59 iso ti«é 172.18 9.92 523.64 
5.610 0.52 15. 67 0.70 171.39 8.64 525.14 
6.934 0.47 15.77 0.66 170. 64 7,56 526.37 
8.715 0.43 16.06 0.60 169.85 6.55 527.51 
12. 497 0. 37 16.79 0.55 168.63 5.16 529.07 
21.923 0.29 1g. 23) ~=«G. a9 166.81 3.34 Geljge 
48. 410 0.19 20.09 0.44 164.74 149" $9045 
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TABLE III-3 
EXACT FLOW SPECIES CONCENTRATIONS FOR RUN 1-CM 


Species Concentration 
(lb-mol/lb) x 104 


Area O H O, Ho OF H 79 
Et) ai ieee 23.80 7203 177.94 PN S18) 507. 86 
1.941 ea 7 Es 24. 00 2.04 L715 20.73 902.15 
2200 Se 23.78 Za00 DG 19.94 510. 10 
ae 74 U0) eo / 23.30 1 eis Ry 1s) Aaah Ss. 98 on L382 
2.414 1.43 Be 3 1.78 WS = She) V7. ae 912.94 
2.711 NES: 21.16 ees. 174.77 Po09 Sia 26 
3.251 Oe 92 L3a7 0 lake 173. 64 13.45 o1S396 
3.631 0.76 es Gao 173.07 es) 7 o2 lea07 
4,055 O02 15.76 0.80 WHZe 39 05 4 i neh 
4.466 Oia 15. 68 Oee7 BAG 7AM Sh Sle 52.3. 66 
9. 431 0.53 ono 0.70 70 8.76 o2io000 1 
6.427 0.48 15.63 0.66 170.90 7.88 326.0: 
JA AS} 0.42 16.08 0.69 169.68 6.30 927. 80 
15.594 0.34 E72 0.52 167.90 4.37 20. ge 
Fateh. S118 0.26 18. 83 0. 47 166.05 2262 53. 910 
46.046 O16 eS, oro L 0.44 164. 87 Way 933. 06 


E6 


34. 


41. 


0.71 


TABLE III-4 


Species Concentration 


(lb-mol/lb) x 104 


67 


162.23 


EXACT FLOW SPECIES CONCENTRATIONS FOR RUN 2-S 


EEE _ —— 


TABLE III-5 
EXACT FLOW SPECIES CONCENTRATIONS FOR RUN 2-C 


Species Concentration 
(lb-mol/lb) x 10 


Area O H O, H, OH H,O 
SEAL ee ee ee eee ae ee 
Ihe, Ihre 3.59 32.14 3.59 183e53 32.48 492.50 
Ils NE 3.08 ch ee as) 3.47 181.88 29.66 496. 06 
lays 29 30805 323 180. 39 27.14 499.35 
ler 22 2.49 28.96 2.97 178. 82 24.54 S020 7 
2. 046 2.18 27.69 27 17ige25 21.94 506.28 
2304 1.91 26.30 2.39 igre 20308 509.03 
2.621 Vas 24. 63 2107 175-11 18.11 511. 86 
3.062 ere) DO OE 1. 66 173.95 lSev2 515.389 
3.445 1.07 20.46 1.38 720 14.00 517 ao2 
4.073 0.93 19.73 lez Wee 32 12.44 519.93 
4.764 0.87 19.79 129k V7 aeo7 11.37 5220 
6.093 0.79 20.10 aos 170. 43 9.86 572093 
9.425 0.69 21.28 0.98 168. 44 7.60 525. 49 
13.424 0.61 22.44 0.93 166. 86 6.00 527. 26 
18.236 0255 23.56 0.90 165.54 4.76 528.64 
29.216 0.45 25102 | 0.85 lGse7 | Bh IY 530. 46 
37.160 0. 40 25096 0.84 162.92 243 Sileegs 
43.092 0.38 26.36 0.83 162. 48 2.06 BoieeS 


TABLE III-6 
EXACT FLOW SPECIES CONCENTRATIONS FOR RUN 3-S 


Species Concentration 
(Ib-mol/lb) x 104 


Area O H O H OH HO 
Al 2 2 2 
1.083 3.80 33.23 3. 80 184.38 33.87 490.48 
1. 212 3. 49 33.24 3.84 183.11 31.80 492.78 
1.444 3.17 32. 42 3. 66 181.16 28.46 496.81 
2.008 2.75 31.94 3.34 178.03 23.82 502.49 
2.625 2.55 32.46 3.19 175.79 20.89 505.93 
3.194 2.44 33.22 3.11 174.22 19.01 508.07 
6.714 2.15 37.58 2.99 168.42 12.85 514.77 
13.303 1.90 42.12 2.95 163. 47 8.13 519.81 
22.980 1.72 45.15 2.95 160. 33 5.24 522.87 
35.282 1.61 46.84 2.96 158.61 3.67 524.53 
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TABLE III-7 
EXACT FLOW SPECIES CONCENTRATIONS FOR RUN 3-CM 


Species Concentration 
(lb-mol/lb) x 104 


Area O H O95 Ho OH H»0 
oo. eee See SL ee eee. 
1.090 3.80 33.23 3. 80 184. 38 33. 87 490. 48 
Yo el 3. 49 33.69 3. 86 182.53 31.00 493.53 
1.401 3.31 33.08 3.78 181.59 2320 495. 67 
1.780 Zao 29.95 3. 16 17s. 90 24.56 902.25 
2.024 2.20 MET ES) Zia Z 177. 48 22. 06 905.95 
2.412 1.69 24.65 2 aeeZ 175.76 18.99 S10. 82 
2.781 1.34 Pa VAY, 70 174.58 Ouray 514. 43 
3.275 Us92 18.53 |i lif 173. 41 13.53 518.96 
3. 676 0.79 17.45 Teno V7. 63 12.18 920.75 
4.023 0.69 16.54 0.89 172.42 11.14 922.14 
4.612 0. 66 EGie7 72 0.85 171. 86 10.39 Seee oO 
5. 406 0.65 17.33 0.85 171.18 9.69 52057 0 
8.911 0.62 19.38 OPE fs) 169.05 7.08 929. 90 
Loeoo2 0.58 Zeal 0.85 GG 7 7 Da AE, 928.18 
Zoe 39 0.595 PAPA e134) 0.85 165.55 4.12 929. 39 
A A AY 0.53 23.44 0.85 164.91 3250 530. 02 
VAS LES, 0.51 23.93 0.85 164. 42 3.04 930.50 
38.576 0.50 24.33 0.85 164.02 2.66 530. 89 
41.708 0.50 24. 49 0.85 163. 86 2.50 Sol 0is 
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Area 
Ratio 


Velocity Pressure Temp 
(ft/sec) (lb/ft?) (Deg R). Number 


| (Ib-sec/lb) (ib/ft3) 


TABLE III-8 


EXACT FLOW VALUES FOR RUN 1-S 


- oe i ¢ RSS —ii - Se FE TOTES TOE, «SS TAIT 


Mach 


3.78 


4.08 


Zi 


Specific 
Impulse 


340.1 


344.2 


351.9 


358.1 


366.1 


372.4 


370s 7 


382.7 


386.8 


391.6 


399.2 


413.3 


420.7 


434.0 


442.9 


449.6 


458.8 


x 103 


12. 06 


11.04 


0.51 


0.33 


Density Gamma 
(Frozen) 


1,864 


1.156 


1,146 


1.148 


i, 154 


NS 


1.163 


1,169 


1.173 


1.178 


1.185 


ie 7 


i 203 


ees 


1.220 


1.226 


1.236 


TABLE III-9 


EXACT FLOW VALUES FOR RUN 1-C 


Area Velocity Pressure Temp Mach Specific Density Gamma 
Ratio (ft/sec) (lb/ft2) (Deg R) Mumber Impulse x 103 (Frozen) 
(lb-sec/lb) (lb/ft3) 
Lag723 8743 7201 5256 Lao 340.0 M05 ela OW 7 
1.965 9197 9803 9088 1.94 348.7 10.05 1.181 
2.234 9589 47 47 4937 2.06 356.4 8.48 1.163 
2.560 9953 3886 4819 2 lig 364.1 Zid 3 lia] 
2.881 10238 3292 4740 2020 370.4 6.16 Ls 
3.253 10512 2786 4666 2.33 376.6 2.31 Lawl? 6 
3.441 10633 2582 4633 2.37 379.4 4.96 1.127 
3.726 10800 2319 4586 2.42 383.2 AS) cali h31 
4.023 10956 2092 4538 2.47 386.9 4.12 1.144 
4.518 11189 1783 4442 290 3921 3.59 . 1.184 
5.610 11599 1321 4242 Zz. 70) 401.3 2.79 1.190 
6.934 11966 988 4051 2.85 409.6 2019 oleae 
8.725 @12627 725 3850 3.01 417.9 1.69 ay 2i02 
12.497 12836 447 3547 3.26 429.7 Jweaeled 1.210 
2. 92.3 13510 212 3107 3.65 445.5 0.61 Lda 
29.235 13809 145 2898 3.86 452.6 OutS  dey220 
48.410 14269 79 2999 4.23 463.4 0.26 1.241 
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TABLE III-10 
EXACT FLOW VALUES FOR RUN 1-CM 


Area Velocity Pressure Temp Mach Specific Density Gamma 
Ratio (ft/sec) (lb/ft2) (Deg R) Number Impulse x10% (Frozen) 


ll el Mb oaeeib) ie 


1.727 8752 7169 5252 1.82 340.2 12.01 1.201 
1.94] gys7 5919) 5100 1.93%  347%9 10.22 W975 
2106 ©9864 5863. 5032 1.99 25128 ama9, dlaba 
2210 986] 401% 4965 2.066 3968640 S266 | Ile 
2.414 9793 4251 4893 2.11 360.9 768 1..4S 
2.711 10087 3602 4808 2.20) 367.3 6.64 1,136 
3.251 10509 2798 4687 2.33 376.7 s4o1 i. wee 
3.631 10746 2408 4618 2.40 382.1 (a55 , els 
4.065 +10973 2074 4546 2.4785 30704 4.08 1.182 
Sn4c6 allg7) 1810. 4463 2.54) 39pm7 3.64 1,185 
6.461 Sins45 1380) 4278 2.6¢m eon 2.90 1.190 
6.427 11843 1096 4118 2.80 406.9 2ase: Aw105 
9,173 12408 676 3805 3.05 419.8 1.60 1.203 

15.594 13122 333 3367 3.42 436.4 0.89 1,215 

28.516 13788 150 2915 3.84 452.1 0.46 1.228 

46.046 14231 0 2587 4,08  4aQ05 0.28 1.240 
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TABLE III-11 
EXACT FLOW VALUES FOR RUN 2-S 


Area Velocity Pressure Temp Mach Specific Density Gamma 
Ratio (ft/sec) (lb/ft?) (Deg R) Number Impulse x10 ; (Frozen) 
(Ib-sec/lb) (lb/ft”) 


em SS SSS SSS SS SS SSS a SSS SSS SSeS SE ew ee Eee, ee 





ies gies l67o7 5650 laasm 3oGno 24.05 1.734 
ie2i7 Meso 13119 5534  ve4em ileus goes7 #1.aB4 
1.481 8190 9141 5338 1.68 32901 14.96 1.138 
e772 Wees5c 6814 5174 1.85 341.8 im s7 wie 
2.198 9525 4891 4980 2.04 355.3 8.67 1.154 
2.417. 9787 4248 4895 2.11 360.8 8.67 1.158 
2.717 10091 3579 4789 2.20 367.2 eac2  lelen 
3.270 10531 2743. 4620 2.34 376.7 527 | 72 
3.636 10765 2360 4523 2.42 381.9 4.64 77 
4.125 11028 1977 4408 2.51 38787 3.99 1.182 
45635 11257 1680 4300 2.60 392.8 3.48 1.186 
Bz46 “aNGS0 124% 4106 2.75 401.6 D.71 12193 
7.160 017 925 3911 2.90 410.0 2.11 1.200 
10.482 12583 553 3587 3.17 423.0 1.38 1.209 
16.422 13154 304 32308 3.46 43aN3 0.84 1.219 
26.911 13688 158 2866 3.83 448.9 0.49 1.230 
34.625 13929 113 2692 4.01 454.6 0.38 1.236 
41.552 14091 89 2571 4.15 458.4 0.31 1.241 
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TABLE III-12 
EXACT FLOW VALUES FOR RUN 2-C 


Area Velocity Pressure Temp Mach Specific Density Gamma 
Ratio (ft/sec) (lb/ft2) (Deg R) Number Impulse x 10% (Frozen) 


ee eS eres el 


lL. ls 6715 15707 9650 1.33 306.0 24.05 1, 580 
NAS 7606 11514 5425 eee 318.6 18.43 VY, 190 
1.473 8195 9102 9272 b.6g 328.6 15.04 1. o/ 
lea722Z 8783 7011 9105 1.84 339.5 12.00 1.160 
2.046 9331 9359 4941 1.30 350.4 9. ou 1.146 
2.314 9674 4468 4846 2.09 357.6 8.11 1.136 
Zi621 9990 3743 4763 7 M8 364.5 6.93 1. 7 
3.062 10353 3017 4670 2e20 872.7 972 1.123 
3.445 10610 2571 4602 Ze 378.6 4.97 1.124 
4.073 10966 2035 4459 | 2.49 386.5 4.06 1.185 
4.764 PLZ 1634 4311 2.60 393.4 3.38 1.189 
6. 093 L725 1163 4086 jay Ti 403.5 2.94 1.196 
9.425 12414 643 3703 3.08 419.2 1.55 1. 207 
13.424 12894 400 3412 3.33 430.4 09 1.214 
18,236 13263 266 3173 od. 54 439.0 0.75 1,221 
29.216 13761 143 2830 3.88 450.7 0.45 1.231 
37.160 13988 104 2664 4.05 456.1 0.35 1. 237 
43.092 14119 85 2565 4.16 459.2 0.30 1.241 


%§ 


Area 
Ratio 


6404 


7262 


8127 


9289 


10020 


10481 


11028 


11855 


12787 


13382 


oe) 7 


Velocity Pressure 
(ft/sec) ' (lb/ft2) 


17277 


13027 


ee RY, 


9415 


3605 


2704 


1853 


954 


379 


183 


103 


TABLE III-13 


Temp 


Mach 


(Deg R) Number 


—— i ee SS SS SS —$—$—$—$——_—_—— SSS ERE AEE O——=——EE = ee ieee 
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EXACT FLOW VALUES FOR RUN 3-5 


Specific Density Gamma 


Impulse 


Si Zul 


312.8 


326.7 


348.6 


363.6 


373.3 


385.3 


403.8 


425.2 


439.0 


448.3 


x 10° 


26.17 
20.61 
15. 46 
Je 3 
G.9 0 
9. 42 
3.95 
AS, 
1.07 
0.59 


0. 37 


(Frozen) 


(lb-sec/lb) (lb/ft3) 


1. 342 


1.169 


PEO 


ee 


V.T8s 


1.194 


1.250 


Area 
Ratio 


11605 


9821 


6631 


9472 


4254 


3486 


2801 


2386 


2108 


1737 


1389 


698 


326 


206 


157 


123 


99 


89 


TABLE III-14 


Velocity Pressure Temp 
(ft/sec) lb/ft2) 


= Yamba = ELS 9 =EQrE ego ee SE Oe 


Mach 


(Deg R) Number 


yy | 


27 


EXACT FLOW VALUES FOR RUN 3-CM 


Specific Density Gamma 


Impulse 


302.1 


317.3 


324.4 


341.1 


349.2 


359.5 


367.3 


375.8 


381.5 


385.8 


392.0 


398.8 


417.6 


435.2 


444.4 


449.4 


453.6 


457.2 


458.8 


x 103 


26.17 
18.71 
16.20 


11. 48 


(Frozen) 


(lb-sec/lb) (lb/ft) 


1.286 
1.168 
1,154 
1.126 
1,120 
Ll. 112 
1. 122 
1 ho) 
1.143 
1.048 
1.203 
1.205 
1,212 
1.221 
1,227 
1.232 
1.236 
1.240 


1.242 
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each area ratio to be chosen to keep the flow near equilibrium. The effectiveness 
of the resulting nozzle protiie is evaluated by comparing exact calculations with 
equilibrium caiculations. Contouring is shown to be effective for small nozzles 
when a correction factor is applied. The hydrogen-oxygen system is used as an 
example, 


FORM 
DD 223. 1473 UNCLASSIFIED 
79 Security Classification 


UNCLASSIFIED 


Security ~ Security Classification = 


Ree MOnDs 






Dissociation 


Optimum; Nozzte, Profiles, Frozen, Flow, 
Expansion, Gas Mixtures, Equilibrium, Recombinatior 








| 
: 


INSTRUCTIONS 


1. ORIGINATING ACTIVITY: Enter the name and address 
of the centractor, subcontractor, grantee, Department of De- 
fense activity or other organization (corporate author) issuing 
the report. 


2a. REPORT SECURITY CLASSIFICATION: Enter the over- 
all security classification of the report. Indicate whether 
‘‘Restricted Data’’ is included. Marking is to be in accord- 
ance with appropriate security regu gop. 


26. GROUP: Automatic downgrading is specified in DoD » 0S la 


rective 5200.10 and Armed Forces Industrial Manual. Enter 

. the group number. Also, when applicable, show that optional 
markings have been used for Group 3 and Group 4‘as author- 
ized. 


2. REPORT TITLE: Enter the complete report title in all 
capital letters. Titles in all cases should be unclassified. 
Uf a meaningful title cannot be selected without classifica- 
tion, show title classification in all capitals in parenthesis 
immediately following the title. 


4. DESCRIPTIVE NOTES: If appropriate, enter the type of 
report, e.g., interim, progress, summary, annual, or final. 
Give the inclusive dates when a specific reporting period is 
covered. 


5. AUTHOR(S): Enter the name(s) of author(s) as shown on 
or inthe report. Enter last name, first name, middle initial. 
lf military, show rank and branch of service. The name of 
the principal author is an absolute minimum requirement. 


6. REPORT DATE: Enter the date of the report as day, 
month, year; or month, year. If more than one date appears 
on the report, use date of publication. 


7a. TOTAL NUMBER OF PAGES: The total page count 
should follow normal pagination procedures, i.e., enter the 
number of pages containing information. 


7b. NUMBER OF REFERENCES: Enter the total number of 
references cited in the report. 


8a. CONTRACT OR GRANT NUMBER: If eppropriate, enter 
the applicable number of the contract or grant under which 
the report was written. 


8b, 8, & 8d. PROJECT NUMBER: Enter the appropriate 
military department identification, such as project number, 
subproject number, system numbers, task number, etc. 


9a. ORIGINATOR’S REPORT NUMBER(S): Enter the offi- 
cial report number by which the document will be identified 
and controlled by the originating activity. This number must 
be unique to this report. 


9b. OTHER REPORT NUMBER(S): If the report has been 
assigned any other report numbers (either by the originator 
or by the sponsor), also enter this number(s). 


10. AVAILABILITY/LIMITATION NOTICES: Enter any lim- 


itations on further dissemination of the report, other than those| 


FORM 


vanes 1473 (BACK) 


DD 





mee by security classification, using standard statements 
such as: 


(1) ‘‘Qualified requesters may obtain copies of this 
report from DDC.’’ 


| 
(2) ‘*Foreign announcement and dissemination of this 
report by DDC is not authorized.”’ | 

| 


(3) ‘*U. S. Government agencies may obtain copies of 
..» this report directly from DDC. Other qualified DDC 
“users shall request thréugh © ~ % ab or dla 


(4) ‘*U. S. military agencies may obtain copies of this 
report directly from DDC. Other qualified users 
shall request through 





(5) ‘All distribution of this report is controlled. Qual- 
ified DDC users shall request through 





99 
b 


If the report has been furnished to the Office of Technical 
Services, Department of Commerce, for sale to the public, indi- 
cate this fact and enter the price, if known 


11, SUPPLEMENTARY NOTES: Use for additional explana- 
tory notes. 


12. SPONSORING MILITARY ACTIVITY: Enter the name of 
the departmental project office or laboratory sponsoring (pay 
ing for) the research and development. Include address. 


13. ABSTRACT: Enter an abstract giving a brief and factual 
summary of the document indicative of the report, even though 
it may also appear elsewhere in the body of the technical re- 
port. If additional space is required, a continuation sheet shall. 
be attached. 


It is highly desirable that the abstract of classified reports 
be unclassified. Each paragraph of the abstract shall end with 
an indication of the military security classification of the in- 
formation in the paragraph, represented as (TS), (S$), (C), or (U). 


There is no limitation on the length of the abstract. 
ever, the suggested length is from 150 to 225 words. 


14. KEY WORDS: Key words are technically meaningful terms 
or short phrases that characterize a report and may be used as 
index. entries for cataloging the report. Key words must be 
selected so that no security classification is required. Identi- 
fiers, such as equipment model designation, trade name, military 
project code name, geographic location, may be used as key 
words but will be followed by an indication of technical con- 
text. The assignment of links, rales, and weights is optional. 





How- 


UNCLASSIFIED 


Security Classification 























